Summary 

A “robust closed-loop system” maintains desired stability and performance when the nominal 
system is subject to plant variations. Multivariable robustness analysis and robust control 
design tools take into account the uncertain nature of the plant models used in the design 
of feedback control systems. A fundamental research issue is to find the means of generating 
descriptions of model uncertainty consistent with robust multivariable analysis and design tools. 
One approach, based on sector stability theory, provides a frequency-domain description of the 
uncertainty associated with a particular model, given bounds on individual parameters in the 
model. This paper explores the application of this sector-based approach to the formulation of 
useful uncertainty descriptions for linear, time-invariant, multivariable systems. A review of basic 
sector properties and the sector-based approach are presented first. The sector-based approach 
is then applied to several general forms of parameter uncertainty to investigate its advantages 
and limitations. The results indicate that the sector uncertainty bound can be used effectively 
to evaluate the impact of parameter uncertainties on the frequency response of the design model. 
Inherent conservatism is a potential limitation of the sector-based approach, especially for highly 
dependent uncertain parameters. In addition, the representation of the system dynamics can 
affect the amount of conservatism reflected in the sector bound. Careful application of the 
method can help to reduce this conservatism, however, and the solution approach has some 
degrees of freedom that may be further exploited to reduce the conservatism. 

1 Introduction 

Uncertainties associated with plant models used in control design are known to limit 
attainable system performance. Attaining higher levels of performance for uncertain multi- 
input-multi-output systems is a key focus area in current controls research. One class of 
multivariable control analysis and design methods explicitly takes into account plant model 
uncertainties. These methods involve the generalization of Nyquist theory to multivariable 
systems by exploiting the properties of matrix norms. For example, methods of accounting 
for unstructured uncertainty use singular values as a basis for robustness tests that assess the 
likelihood of a system remaining stable in the face of particular uncertainty representations. 
These tests, however, require knowledge of the uncertainty magnitude over all frequencies to 
assess robustness accurately. 

Consider the graphical interpretation of a typical robustness test for a linear system with 
unstructured multiplicative uncertainty reflected at the plant output depicted in figure 1 
(Postlethwaite and Foo 1984). The upper curve is the reciprocal of the maximum singular 
value of the complementary sensitivity function. The lower curve is the maximum singular value 
of the uncertainty. The condition for robust stability is that the two curves do not intersect. The 
smallest distance between the two curves determines the level of robustness (i.e,, f the robustness 
margin). The frequency-dependent characterization of the uncertainty is therefore critical in 
accurately assessing system robustness. Other methods, such as those based on structured 
singular values, use the knowledge of the uncertainties in the form of weighting functions for the 
same reasons. 

Attempts to use the complementary sensitivity function alone as a test for robustness are 
likely to give erroneous and misleading results. This is the case in figure 1. Point a corresponds 
to the point of minimum robustness when the uncertainty is assumed to be independent of 
frequency. Point b, on the other hand, corresponds to the point of minimum robustness when the 
frequency dependence is considered. A true test for robustness requires both curves. The upper 
curve alone only provides a relative robustness measure at best and an optimistic indication of 
the level of robustness and an erroneous indication of the critical frequency range at worst. (A 
nomenclature section appears after the references.) 
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Figure 1. Robustness test for unstructured uncertainty. 


Design methods that explicitly take into account plant model uncertainties generally assume 
that descriptions of the uncertainties, for example, in the form of additive or multiplicative 
transfer function matrix blocks, are available as part of the plant model (Doyle and Stein 
1981; Lehtomaki et al. 1984; Postlethwaite and Foo 1984; Ridgely and Banda 1986). Such 
representations of plant uncertainties are not usually part of the design model and methods of 
representing them are not well developed (Anderson and Schmidt 1989; Stein and Doyle 1978; 
Kwakernaak 1985). A fundamental research issue is finding the means to generate descriptions 
of plant model uncertainty consistent with robust multivariable analysis and design tools. 

One approach to formulating uncertainty descriptions has been developed by Safonov (1978) 
and is based on sector stability theory (Zames 1966a, 1966b). The sector-based approach is very 
general and permits many forms of uncertainty (both linear and nonlinear) to be considered. 
It is also completely compatible with multivariable systems and singular-value-based robustness 
measures. The fundamental goal of the approach is to determine a bound on a given function 
of parameters for which values are within known bounds. The function may, for example, 
correspond to the frequency response of the plant, and the bounded parameters may correspond 
to various plant model quantities that, while uncertain, are known well enough to place a bound 
on their admissible values. 

There are, however, many aspects of the sector-based approach that need to be better 
understood before it can be used with confidence. This paper employs sector stability theory 
concepts using Safonov’s (1978) approach to develop an uncertainty description in a context that 
is consistent with linear, time-invariant, multi-input-multi-output (LTI-MIMO) systems. This 
study identifies key properties of the sector-based approach and highlights important implications 
of its application to robust control analysis methods. It is important to note that although 
this paper reviews the theoretical development performed by Safonov and applies it to linear 
multivariable systems, the emphasis is on the application (i.e . , the specific issues associated with 
the application process and the resulting implications). 

The paper begins with a review of the sector concept and the basic sector properties 
required to formulate the sector-based uncertainty representation in section 2. Section 3 
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begins by presenting a block diagram representation of uncertainty for LTI-MIMO systems 
that is consistent with the sector theory. The sector theory is then combined with the 
uncertainty representation to result in a sector-based uncertainty description that is compatible 
with multivariable robustness analysis and design methods. In section 4, this compatibility is 
demonstrated by combining the sector-based uncertainty description with one of the singular- 
value-based robustness conditions (Doyle and Stein 1981; Postlethwaite and Foo 1984). Finally, 
in section 5, the sector-based approach is applied to models that exhibit three classes of parameter 
uncertainty. The first class of uncertainty is associated with a MIMO plant model in which 
only parameters that appear in the numerators of the system transfer function matrix are 
uncertain. The second class of uncertainty is associated with a MIMO plant model in which 
only parameters appearing in the transfer function matrix denominator are uncertain. These 
examples illustrate a method for computing the multivariable sector bound and help identify 
and emphasize some key properties and limitations of sector-based uncertainty descriptions. 
The third class of uncertainty is for a single-input-single-output (SISO) plant model that is 
representative of real-world uncertainties that occur in the transfer function numerator and 
denominator simultaneously. 

2 Sectors and Their Properties 

The term sector comes from a geometric interpretation of some basic concepts of the stability 
theory from which sectors were developed (Zames 1966a, 1966b). These concepts are reviewed 
to understand the value of sector theory in describing plant model uncertainty and to identify 
its limitations. 

Consider a real-valued function / [u(t), C] of the real variable u(t) and parameter the graph 
of this function is shown in figure 2. The variable u is itself a function of the independent variable 
t, which might typically be time. The graph of the function lies between the dotted lines with 
slopes labeled a and b. The region of the plane between the dotted lines containing the graph 
of the function forms a sector of the plane. If the function is known to lie inside the sector for 
all values of the parameter then the sector is a bound on all values of / [u(t), £]. As a result, 
/[u(/),£], as a function of u(t), can be characterized in terms of an expression involving the 
slope of the lines bounding the sector. 


Figure 2. Geometric interpretation of sector bounds. 
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to be the “center” and the “radius” of the sector, respectively. The fundamental property of 
the relationship between the function /[u(/),£] and the sector in which it is contained can be 
described mathematically in terms of c and r as follows: 


!/[«(<),(] -cu{t) I < |ru(/) I (2) 

In addition, equation (2) can be integrated over the range of / to give a sector bound involving 
the 2-norm, or 

11/ [«(<),(] ~cu(t) || 2 < ||ru(/) ||2 (3) 

where the 2-norm for a real-valued scalar function is defined by 

11/ (/) ib — r fa) fa) dt ( 4 ) 

j — 00 

Any function that satisfies equation (3) is said to “lie in the interior of the sector with center c 
and radius r.” 

The concept of the sector can be extended to the complex domain (s-plane) by applying 
Parseval’s formula (Papoulis 1962). ParsevaPs formula equates the inner product of a real- 
valued function to the inner product of its Fourier transform. That is, 

/ OO 1 TOO 

/ it) f it) dd — — /l rl/lrl (U (5) 

-CO ^ J — CO 


where / is the Fourier transform of f, 

/ OO 

exp {—jut) f (t) dt (6) 

-OO 

Therefore, in the complex domain the sector bound is simply 

11/ (s, () - C u (s) || 2 < 1 1 rO. (s) || 2 (s - ju>) (7) 

where the 2-norm for a complex-valued scalar function is defined by 

Il/Hll2= r ./"'-■ )/!.• I doj (8) 

J — 00 

Now consider a mathematical generalization of the concept of a sector. Allow f(s) to 
be F(s)u(s), where u(s) is a complex-valued vector and F(s) is a matrix of complex-valued 
functions. Then allow c and r to be expressed as matrices of complex-valued functions C(s) and 
R(s), respectively. The generalized form of the sector bound equation (eq. (7)) is then written 
as 


[F(s)-C(s)]u(s)|| 2 < ||R(s)u(s) || 2 (s = /w) (9) 
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The physical interpretation of the sector is now obscured, but the same terminology is used 
when discussing the properties of equation (9), that is, center, radius, and sector bound. The 
sector bound in equation (9) is the form that leads to the application of sector theory to 
uncertainty modeling. However, some additional manipulations are required to relate the sector 
bound to recently established robustness concepts. The properties required to accomplish this 
are presented below. 

2.1 Property 1 

A significant advantage of the sector-based approach to bounding matrices of complex- valued 
functions is that the bound can be related to singular-value properties of the matrices. The 
2-norm for functions of the type A(s)x(s), where A(s) is a complex matrix and x(s) is a 
compatible vector, has the form 

|| A (s)x (s) || 2 = — / x*(s)A*(s)A(s)x(s) ds (10) 

where the integrand is the Euclidean norm of A(s)x(s), denoted by ||A(s)x(s)||_g;. By writing 
equation (9) in the form of equation (10) it can be readily shown that if 

1 1 F (s)u(s) - C (s)u(s) || E < ||R'(s.)u(s) \\e 0 = ju>) (11) 

for all 1 0 G [0, 00 ), then equation (9) is satisfied as well. That is, if equation (11) is satisfied 

for every value of the independent variable s, then the integral of the left side is less than the 

integral of the right side, This is exactly the requirement for equation (9) to be satisfied. 

Note that the minimum and maximum singular values of a matrix are related to the Euclidean 
norm by 

(12a) 

<L(A)=$'„11^A (12b) 

\\ x \\\e 

where <j(A) and <j(A) indicate the maximum and minimum singular values of the matrix A, 
respectively (Ridgely and Banda 1986). Applying this property to equation (11), one can see 
that equation (9) is guaranteed to be satisfied if 

a[F (s) - C (s)] < a[R(s)] (s - jui) (13) 

for every ui G [0, 00 ) and ||x 0)11 E ^ ; 0. That is, if an upper bound on the left side of 

equation (11) is less than a lower bound on the right side, then equation (11) is satisfied. This 
form of the sector bound can be readily combined with established robustness properties, as 
shown subsequently. 

2.2 Property 2 

Consider the sector bound 


11/ [«(<),(] - cu $) II 2 < IMO II 2 (14) 

If f[u(t), C] can be represented by [c + 6c(t)]u(t), where 6 c(t) represents a perturbation about, 
c, then equation (14) simplifies to 


\\6c{t) u{t) || 2 < 1 1 ru (/ ) || 2 


(15) 
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One interpretation of this result is that if c + 6 c(t) lies in a sector with center c and radius r ? - 
then 6 c(t) lies in a sector with center at the origin (zero) and radius r. This indicates that the 
sector bound can be represented in a form that is independent of the nominal values of the 
parameters. In terms of the graphical interpretation of the sector (as in fig. 2), this corresponds 
to rotating the sector about the origin until it is symmetric with respect to the independent 
variable, as shown in figure 3. The sector bound for c + 6 c(t), equation (14) and figure 3(a), is 
therefore completely equivalent to the sector bound for 6c(t), equation (15) and figure 3(b). 


(a) Sector for c + 6c(f). 


(b) Sector for 8c{t). 


Figure 3. Equivalence of sector bounds. 


2.3 Property 3 

Another key property of sectors is that if there is a collection of sectors that are independent, 
they can be combined into a single matrix sector by concatenating them diagonally. Consider 
the collection of perturbations c + <5cj, which lie in sectors with centers C; and radii iq such that 
none of the c 8 - + <5cj depend on the other terms cj and 6cj (j ^ i). Safonov and Athans (1978) 
have shown that, in this case, 


(C + <5C)x-Cx|| 2 < || Rx || 2 


(16) 
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where 


C = diag [ci, c 2 , c 3 , ...] 

(17a) 

(5C = diag [<5ci, 6c 2 , <5c 3 , ...] 

(17b) 

R = diag [?q, r 2 , r 3 , ...] 

(17c) 


and x is a compatible vector. This fact is important because it allows sector bounds on individual 
system parameters to be used to bound the entire system. 

3 Sector Bounds for the System 

The sector properties discussed above can be combined with an appropriately developed 
system representation to formulate a sector bound for the plant when the variations in its 
parameters have known sector bounds. The following development of the bound parallels the 
proof of Safonov’s sector bound theorem (Safonov 1978; Safonov and Athans 1978, 1981). It is 
less rigorous than the proof, however, and is only intended to explain the key concepts. 

Assume that the plant model is given such that 

y ( s ) — T (s) u (s) (18 ) 

where y(s) and u(s) are the system outputs and inputs, respectively, and T(s) is a linear, 
time-invariant s multivariable operator with uncertain parameters. The uncertainties cause 
perturbations in the plant response relative to the response of the nominal plant model. The 
response of the plant can therefore be written 

y (s) - L yu (s) u (s) + d (s) (19 ) 

where L 2/ . u (s) is the transfer matrix for the nominal plant model and d(s) is the perturbation 
response. The perturbation response is directly dependent on the variations in the plant 
parameters due to the modeling uncertainty. The functional relationship between these 
quantities can be expressed as 

d (s) = f (s, C) u (s) (20 ) 

where f(s, () represents perturbations in the system dynamics due to the parameter variations 
associated with uncertainties in the nominal plant model. 

Assume that each uncertain parameter (j is independent and lies in a sector centered at 
c i with radius iq. Then according to property 3, the uncertainties can be combined into a 
perturbation matrix C + <5C, with the sector bound given in equation (16). By property 2, an 
equivalent sector bound is 

||6C x|| 2 < || Rx || 2 (21) 

In order to exploit this sector bound, the perturbation matrix <5C must be isolated from the 
dynamics. That is, d(s) must be able to be expressed in the following form: 

d (s) = 6L yu (s, 6C) u (s) (22) 

where <5C appears explicitly. The solution is to develop a feedback structure around <5C that 
admits all forms of parameter variations associated with LTI-MIMO systems. 

A structure that allows <5C to be isolated from the plant dynamics is shown in figure 4. 
The perturbation matrix <5C appears explicitly and the dynamic effects are characterized by 
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the transfer function matrices L eM (s), L e „(s), and L 2/ „(s). This structure is similar to that 
developed by Safonov (1978). The elements in the dashed box represent the portion of the system 
response associated with plant uncertainties and are referred to as the uncertainty dynamics. 
The transfer matrix L eM (s) transforms the vector of system inputs u(s) into one compatible 
with <5C. The matrix L 2/ „(s) transforms the output of the perturbation block v(s) into the 
perturbation response d(s). The matrix L e „(s) is required to allow the parameters in <5C to 
appear in the denominator of the system transfer matrices. Without L e „(s) it would only be 
possible to characterize the effects of uncertainties appearing in the numerator of the transfer 
matrix. This fact becomes evident when one attempts to characterize a denominator uncertainty 
without using L e „(s), an exercise that is left for the reader. 


Figure 4. General structure for uncertainty dynamics. 

Consider the block diagram of the system representation in figure 4. It is evident from the 
diagram that 

y ( s ) — ^yu («) u (s) + h yv (s) v (s) (23a) 

e(s) = L e „(s)u(s) + L e „(s)v(s) (23b) 

where e(s) and v(s) are the input and output of 6 C, respectively. Since the perturbations are 
assumed to lie in a sector centered at the origin with radius R (i.e., eq. (21) with x replaced by 
e(s)), the sector bound for <5C can be written 

I|v(s) \\e < ||R(s)e(s) \\e (24) 

The objective is to obtain a sector bound on the overall system. This can be accomplished 
by replacing v(s) and e(s) in equation (24) with expressions in terms of u(s). The remainder 
of this section is devoted to the mathematical development required to isolate u(s). The main 
mathematical tool is the inner product defined by 

(xi(s) ,x 2 (s)) = xf(s)x 2 (s) (25) 

where xj(s) and x 2 (s) are compatible complex-valued vectors. The following discussion is limited 
to operations in the frequency domain and, for the sake of clarity, the argument s is deleted. 

Substituting equation (23b) into equation (24) and expressing the Euclidean norm in terms 
of the inner product IK*) III = {(•) , (•))) results in 

(v , v } ^ (L ew u + L ev v) , R (L ew u + L ev v )) (26) 
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which can be expanded with inner product operations to 


0 > (v, (I - L*„R*RL e „) v) - (u, L* M R*RL e „v) 

- (v , L*„R*RL e „u) - (u, L*„R*RL e „u) (27) 

Define the following: 

Q„ = I - L*„R*RL e „ (28a) 

R„ = L* u R*RL eu (28b) 

S v = L*„R*RL e „ (28c) 

Substituting these definitions into equation (27) results in the simpler expression 

0 > (v, Q„v) - (S„u, v) - (v, S„u) - (u, R„u) (29) 

Now define 

T® = Q)T 1 S„ (30) 

and note that Q„ usually has an inverse 1 and that Q* = Q„. Also note that the definition of 
T„ has the following properties: 


<(v -T„u),Q„(v -T„u)) = (v, Q„v) - (u, T*Q„v) 

- (v, Q„T„u) + (u, T*Q„T„u) (31) 

and 

T*Q„ = = S* (32) 

Therefore, substituting T„ into equation (29) and using the properties in equations (31) and (32) 
allows equation (29) to be simplified to 

((v-T„u),Q 0 (v-T„u)) < (u, P„u) (33) 

where 

P„ = R„ + SJQ^S,, (34) 

In order to eliminate v from equation (33) a special property of matrices is employed. Let A 

and B be appropriately dimensioned complex matrices, with A being positive definite. If A is 
defined as 

A = (BA-'Bj " 1 (35) 

then 

(x,B*ABx) < (x,Ax) (36) 

This is lemma A1 (Safonov 1978) and it is proved in appendix A. 

Let A = Q„ and B = h yv ] then, 

(x, L yv (liy V Q v Lj/^x) < (x, Q v x) (37) 

1 A condition for the existence of Qp 1 is discussed in section 5.1. 
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If x in equation (37) is replaced by (v — T„u) and compared with equation (33), then one finds 
that 

(L y V (v — Tj,u) , ^-‘yv^j Ly V ( v — T^u)) < (u, P J,u) (38) 

Substituting equation (23a) into equation (38) results in 


( [y (kj/M + Ly V Tj,) uj , ^yv^j [y i^yu d~ Ly V Tjj) u] ) < (u, P„u) (39) 

Note that equation (39) is satisfied whenever equation (24) is satisfied. All dependence on e 
and v has been removed. The information about the sector bounds on the system parameters 
is now completely contained in P„, Q„, and T„. All that remains is to simplify equation (39) 
into a form consistent with the sector bound equation. 

Define 

T n — L yu + L^T*, (40a) 

Qra = 1 (40b) 

Pra = Pv (40c) 

Substitution of equations (40) into equation (39) results in 


((y - T„u) , Q„ (y - T„u)) < (u, P„u) 


(41) 


Now, since Q ra and P ra are Hermitian matrices, they can be factored 2 so that equation (41) can 
be written as 

||Q™ /2 (y-T ra u)|| 2 < ||Py 2 u|| 2 (42) 

Substituting equation (18) into equation (42) results in 

||Qw / 2 (T — T ra )u|| 2 < 1 1 P i^ /2 u| 1 2 (43) 


where, in summary, 


T n — L yu + L yv (I — L*„R*RL e „) 1 L*„R*RL e 


P„ = L* M R* (I - RLggL^RT 1 RL 


Q? 


L V v (I - L*„R*RL e „) _1 L* 


eu 

1-1 


yv 


(44a) 

(44b) 

(44c) 


and where the expression for P ra was obtained by use of the matrix inversion lemma (pp. 48-49 
of Brogan 1974). 

Equation (43) is in the form of a sector bound that is slightly more general than that discussed 

1/2 

previously. (See eq. (9).) Notice that the left side has the added weighting Q ra ' , which did 
not appear in the earlier form. This causes little difficulty, as shown subsequently. Therefore, 
equation (43) is a generalized sector bound on the overall plant response given the sector bounds 
on the uncertain plant parameters. 

2 The factoring method used in the examples is the Cholcsky decomposition(Strang!980). 
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4 Robustness Analysis Using Sector Bounds 

The sector bound shown in equation (43) can be used to compute an upper bound on T 
(i.e. , an uncertainty bound) that is consistent with LTI-MIMO robustness analysis and robust 
control design methods. Begin by considering a singular-value form of the sector bound from 
equation (43): 


<7 


Qra /2 (T 



< <7 



(s-ju-,u G [0, oo)) (45) 


If equation (45) is satisfied, then equation (43) is satisfied as well. Equation (45) can be expanded 
through the use of singular-value properties (Ridgely and Banda 1986) as follows: 


a 


Qra /2 (T 


T„) 


< 


1/21 

n 


pl/2 

1 n 


<= cr [ T - T„] < a 

since, for compatible arbitrary matrices A and B, 

<j [AB] < a[A]a[B 


q ,; 1/2 1 


(46) 


and for dehnite A, 


-1 


1 

^[Aj 


Thus, if equation (46) is satisfied for all s — jui (u G [0, oo)), then the sector bound (eq. (43)) 
is also satisfied. Equation (46) is a general singular-value representation of a sector-based plant 
model uncertainty bound. It is in a form that can be readily combined with recently established 
singular- value-based robustness tests. 


The combination of the uncertainty bound with a robustness test requires the choice of a 
specific characterization of the plant model uncertainty. For the purposes of demonstration, 
consider a typical formulation with additive uncertainty as shown in figure 5. This uncertainty 
formulation is used throughout the remainder of this paper. Other uncertainty formulations 
(e.g., multiplicative uncertainty at the plant input or output) give similar and analogous results. 


Figure 5. Additive plant uncertainty. 
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A condition for guaranteed stability of the closed-loop system depicted in figure 5 subject to 
the uncertainty A is 

a [A] < — f — r {s-jiO]LO G [0, oo)) (47) 

W [I + KG] -1 K 

where G is the nominal plant transfer matrix and K is the feedback compensation (Doyle and 
Stein 1981; Barrett 1980; Postlethwaite and Foo 1984). The plant uncertainty represented by A 
may not actually be known, but a quantitative upper bound on <j[A] can be obtained with the 
sector-based approach if the uncertain parameters in the plant model have known bounds. The 
objective is therefore to obtain an upper bound on <j[A] from the sector bound equation. 

Simple block diagram manipulations indicate that the open-loop system representation can 
be written as 

T = G + A (48) 

A similar expression for T ra in equation (46) can be written: 

T„ = G+A (49) 

where G corresponds to the nominal plant transfer matrix, L yu in equation (44a), and A 
corresponds to the remaining term in equation (44a). Note that A is a purely mathematical 
quantity that in this case is analogous to A. Substituting these expressions into equation (46) 
results in the equivalent sector bound equation 

<t A — A < <j P ^ 2 <j Q ra 1//2 (50) 

Now consider the single-input-single-output form of equation (50) where A and A are complex 
scalars and P ra and Q ra are replaced by the scalars p and q. Equation (50) takes on the simplified 
form 

|A - A| < | (p/g) 1/2 | (51) 

A geometric interpretation of this equation indicates that A lies inside a circle centered at A 
with radius | ( pjq ) 1//2 |, as depicted in figure 6. 


Figure 6. Geometric interpretation of sector bound. 
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Let a real parameter 6 be chosen such that 


<5 = |A| + | (p/q) 1 / 2 | (52) 

Clearly 6 is an upper bound on |A| (1. e. , 6 > |A|). Generalizing this to matrices implies that if 

6 — a A + <t P^ 2 <7 Q„ 1//2 (53) 

then 6 is an upper bound on <j[A], which can be proved as follows. Let a be a matrix such that 
<j [a] = 6. Then, 

a [a] — a A —a P^ 2 <7 Q ra 1//2 (54) 

=y g_ a - A > a P l/ 2 a Q „ 1 ^ 2 (55) 

since for compatible arbitrary matrices A and B (Ridgely and Banda 1986), 

<t[A - B] > a [A] - a[B] 

Combining equation (55) with equation (50) gives 

<j a — A > a P^ 2 <7 Q ra 1//2 > a P^ 2 g_ Q ra 1//2 > a A — A (56) 

As a result, a is a bound on the admissible uncertainty. This implies that 

<7 [A] < <7 [a] = <5 (57) 


<7 [A] < <7 A + (7 P^ 2 (7 Qra 1//2 (58) 

The right side of equation (58) is a sector-based uncertainty bound that is compatible with the 
singular-value-based robustness condition. The same expression results for other uncertainty 
representations (e.g., multiplicative uncertainty reflected at the plant input). The only difference 
is that A has a different interpretation that depends on the uncertainty representation. 

Equation (58) can be used in the robustness condition, equation (47), in place of cf[A] as 
shown in equation (59) to provide a robustness test that can be readily computed : 


<7 A + (7 P^ 2 (7 Q ra 1//2 < 


cF [1+ GKl -1 K 


Equation (59) can be used as a sufficient robustness condition to verify that known sector 
bounded uncertainties will not destabilize the closed-loop system. The implications of equa- 
tion (59) are the same as those of equation (47). That is, if equation (59) is satisfied, then 
the closed-loop system is guaranteed to be stable for the modeled plant variations. Violation 
of equations (47) and (59) does not, however, imply that the closed-loop system is necessarily 
unstable (Doyle and Stein 1981). 

An interpretation of this result is illustrated in figure 7. The upper curve is a plot of the 
reciprocal of the maximum singular value of [1+ GK] -1 K. The lower curve is a plot of the 
maximum singular value of the “exact” uncertainty A. The middle curve is a plot of the 
right side of equation (58), and since it is an upper bound on the actual uncertainty, it is 
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used in equation (59) to test for system stability subject to the sector bounded uncertainties. 
Equation (59) is satisfied as long as the upper and middle curves do not intersect. This condition 
is potentially more conservative than equation (47), however, since equation (59) may be violated 
while equation (47) is still satisfied. In graphical terms, the upper and middle curves may 
intersect without causing the lower and upper curves to intersect. 


Figure 7. Interpretation of sector-based robustness condition. 


5 Limitations and Applications 

In this section limitations in using the uncertainty bound in equation (58) are discussed. Some 
of the limitations are directly observable from the development itself while others become more 
apparent when the approach is applied to typical forms of parameter uncertainties. Example 
problems are, therefore, used to illustrate these limitations and at the same time to demonstrate 
the solution technique. 

5.1 Limitations 

A couple of conditions on the applicability of the uncertainty bound are linked to equa- 
tions (44). In order for the system sector bound to exist the inverse of 

Q„ = I - L*„R*RL et , 

(eq. (28a)) must exist since T ra ,P ra , and Q ra from equations (44) all depend on the existence 
of Q” 1 , either directly or indirectly. In addition, Q ra only exists if L^Q^L*,, has an inverse. 
Simple tests for the existence of these inverses are the following: 

1. Q„ will have an inverse if <r[RL e „] < 1. (Note that this is only a sufficient condition.) 

2. L^Q^L*,, will have an inverse if Q„ has an inverse and if L yv is an to X n matrix, where 
to < n (to is the number of outputs and n is the number of uncertain parameters) and L yv 
is full rank (i.e. , rank to). 

Proofs of these conditions are presented in appendix A. The consequences of these conditions are 
fairly minor. A violation of the first condition implies that sector radii are sufficiently large that 
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the maximum singular value of the matrix product RL e „ exceeds unity. This can be remedied 
by using a different set of physical units. For example, assume a parameter is represented in 
units of degrees per second and the sector radius is 100. Converting the units to radians per 
second reduces the sector radius to 0.556, which is likely to cause the existence condition to be 
satisfied. 

The second condition is violated when there are less uncertainties than outputs or when 
an output of the system is not affected by any of the modeled uncertainties. It is not likely 
for realistic systems to contain less uncertainties than outputs, since there are typically a large 
number of uncertain parameters. However, if there are fewer uncertain parameters than outputs, 
the parameter must be repeated until the condition is satisfied. In the latter situation, if a 
transfer function in the system does not have any uncertainty, it should not be included in the 
computation of the sector bound. Therefore, the remedy is to reformulate the problem such that 
only outputs that are dependent on uncertain parameters are considered. 

A slightly stronger limitation of the sector theory is that it only applies when the true plant is 
stable for all perturbations inside the parameter sectors. This is because the relationship between 
the frequency-domain and the time-domain representation of the system (via Parseval’s formula) 
is only valid for stable systems. This is not too disconcerting for robustness applications because 
the robustness conditions (e.g., eq. (47)) effectively limit the perturbations to be nondestabilizing 
(open loop) since the conditions require that the number of right half-plane poles of the open-loop 
system not change under admissible perturbations (Doyle and Stein 1981). 

Yet another limitation is associated with the conservatism introduced by sector bounds. In 
the frequency domain, if a single parameter a is known to lie in a sector with center f3 and 
radius p, then a can take on any complex value inside the circle, as shown in figure 8. Thus, 
even if a is known to be a strictly real-valued parameter, bounding it within a frequency-domain 
sector admits complex-valued perturbations 6a. As a result, conservatism is introduced into the 
system sector bound when uncertain parameters are strictly real-valued parameters (or otherwise 
conservatively bounded within a circle). Additional conservatism is introduced by the singular- 
value form of the sector bound since it is obtained through use of inequality properties of singular 
values. 


Figure 8. Sector bound for real parameter. 


Conservatism is also introduced because of the restrictions on cross coupling of the uncertain 
parameters. Recall that by property 3 the uncertain parameter matrix 6C is required to be 
diagonal. This allows the individual parameter sectors to be combined into a single-matrix 
sector. (See eqs. (16) and (17).) However, this also requires that each uncertain parameter 
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be treated as if it is independent of the other uncertain parameters. Also, the fact that <5C is 
diagonal limits the ways by which the system can be represented to formulate the system sector 
bound. This conservatism can best be illustrated by examples. 

The following examples serve to identify some of the factors that influence the conservatism of 
the sector-based uncertainty bound. In addition, the examples demonstrate a method of solving 
for the uncertainty dynamics matrices for general plant parameter uncertainties. 


5.2 Numerator Uncertainty for a Multivariable System 


In this example a simple two-input-two-output system is considered to have the transfer 
matrix 


G (s) — 


1 


■s' 2 + ds + e 


s + a 
s + b 


s 2 + as + c 


(60) 


where *ach of the numerator coefficients is uncertain but is known to lie in a sector. That is, 


a = a+a' b = b + b' c = c + c' (61) 

where, for example, a is the nominal value of a and a 1 is the perturbation of a. The sector 
bounds on the coefficients are such that the perturbation value is centered at the origin with 
radius r (e.g., a 1 is in a sector centered at zero with radius r a ). 

Denote the system inputs by u(s) and the system outputs by y(s). The system frequency- 
domain input-output equation can then be written as 


s“ + ds + e ) y (s) = 


s + a a 

s + b s 2 + as + c 


u (s) + 


a' a' 
b' a's + c' 


u (s) 


(62) 


The perturbation response (as in eq. (19)) due to uncertainties is simply 


d (s) = 


1 


s 2 + ds + e 


a' a' 
b' a's + c' 


u (s) 


(63) 


Block diagram manipulations of the uncertainty dynamics, 6 L yu in figure 4, indicate that the 
disturbance can also be expressed as 


d (s) = L yv (s) SC L eu (s) u (s) (64) 

Notice that L e ;,(s) is zero because there are no denominator uncertainties. 

Combining equations (64) and (65) gives the equation that governs how <5C , L e . Ui (s), and 
Ly(,(s) may be chosen: 


L yp (s) 6C L, , t (s) — 


1 


J 


J 


s 2 + ds + e 


a a 
b' a's + c' 


(65) 


The solution of equation (65) is not unique. There is considerable freedom in choosing the 
transfer matrices in the uncertainty dynamics since the scalar equations that govern the choice 
of the matrix elements are underspecified. The governing scalar equations are also nonlinear, a 
fact that complicates their solution. 

The first step in solving equation (65) is the choice of the <5 C matrix. The simplest way of 
choosing <5C is to assume that each appearance of an uncertain parameter is associated with a 
different sector. For example, since the parameter a appears three times in the input-output. 
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equation, each of the three is assumed to be a different parameter and so is associated with 
an independent sector. The perturbation matrix and the associated sector radius matrix are 
therefore 


<5C = diag [V, a' , a' , b' , c/] (66a) 

R = diag [r a , r a , r a , r b , r c ] (66b ) 


One would expect the system sector bound resulting from this choice for 5C to be overly 
conservative because, in effect, a is allowed to take on a different value every time it appears in 
the system input-output equation. This issue is addressed subsequently. For now, consider the 
process by which L eUi (s) and L^s) can be chosen subject to the choice of <5C. Let 


L^.;, (s) 


Nj/r (s) 

dyp (s) 


Le« (s) — 


(s) 
deu (*) 


(67) 


where N yv (s) and N eu (s) are matrices of numerator polynomials and d yv (s) and d eu (s) are 
denominator polynomials. Recall from equation (23a) that 


y (s) = L yu (s ) u ( s) + L yv (s) v ( s) 

Thus, if dy V (s) is chosen to be the nominal system denominator polynomial, then both L 2/ . u (s) 
and L yv (s) will have a common denominator. Therefore, let 

dyp (s) — + ds + 6 (68) 

Equation (63) and the choice for d yv (s) then require that 


deu (*) — 1 


(69) 


For dimensional consistency N yv (s) must be a 2 X 5 matrix of polynomials and N e . u (s) must 
be a 5 X 2 matrix of polynomials. Therefore, let 


N yv (s) 


h-2 hs h4 h5 

^21 h2 ^23 ^24 ^25 


(70a) 


(®) 


m n m- 2 i m 31 m 41 m 51 
m 12 m 2 2 m 32 m 42 m 52 


(70b) 


Substituting equations (70) and (66a) into equation (65), expanding the elements, and grouping 
like terms results in a system of nonlinear equations for each uncertain parameter. This is true 
regardless of how many times a parameter appears in the system equations. Therefore, there are 
three systems of equations to solve in this example, one associated with each of the uncertain 
parameters a, b, and c. 


System for a: 

1 = l n m n + /i 2 m 2 i + ^i3 m 3i ' 

1 = /nm 12 + l\2 m 22 + ^13 m 32 

0 = ( 2 im n + / 22 m 21 + / 23 m 31 

S = ^21 m 12 + ^22 m 22 + ^23m 32 . 


(71a) 
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System for b: 

0 = (i4?T!-4i 

04 

II 

0 


1 = ^24 m 41 

04 

04 

■»»» 

II 

O 

System for c: 

0 = / 15 J 77.51 

0 — / 15 J 77.52 


0 = l-25 m 51 

1 = /25 m 52 


(71b) 


(71c) 


The solutions to these systems of equations are not unique. However, since the elements of 
N y V (s) and N eUi (s) are numerator polynomials, practical restrictions can be put on them that 
make the solutions tractable. Require the elements m ^ and n ij to be minimal polynomials (i.e., 
they only involve either nonnegative integer powers of s or zero) and to be of the lowest possible 
order and least possible number of terms. For example, m jj or n ij will not be allowed to involve 
terms like s -1 or s 0 ' 5 . 


With these restrictions, reasonable solutions of the systems of equations for nonrepeated 
parameters, such as b and c, are quite simple. A solution to the system for b is 


^14 = 0 ^24 — 1 m 4i = 1 J71-42 = 0 (72) 

and a solution to the system for c is 

^15 = 0 ^25 — 1 m 51 = 0 m 52 = 1 (73) 

The solution to the system of equations for repeated parameters, such as a, is slightly more 
involved since, even with the restrictions, there are many more possible solutions. The system 
for a can be solved by making some “inspired guesses.” Let 

(ll = ^22 = 1 (l2 = hi — 0 (l3 = ^23 = 0 m 22 = s (74) 

The remaining terms in the system for a then become 


mn — m l 2 — 1 m 21 — 0 m 31 — m .32 = Arbitrary 


(75) 


Since m 3 i and 77732 are arbitrary let them both be 1, for example. With these values for the 
numerator terms N yv (s) and N eu (s), the complete uncertainty dynamics representation becomes 


6C — diag [a', a', a', b', c' 



10 110 

1 s 1 0 1 


(76a) 


L</;, (S) 


1 

s 2 + ds + e 


1 

0 


0 0 0 0 

10 11 


L e< , (s) — [0]5x5 


(76b) 


with the associated sector radius matrix R from equation (66b). 

Note that the third column of L^s) is zero, and so the third appearance of the perturbation 
a 1 does not influence the plant responses. This implies that the third appearance of the 
perturbation a 1 is superfluous and may add unnecessary conservatism to the resulting system 
sector bound. It can be dropped from <5C without affecting the ability to appropriately represent 
the uncertainty dynamics tfL^s). 

A key question associated with the solution to equations (71a) presented above needs to be 
addressed: Is there actually additional conservatism in the sector bound associated with an 
excessively repeated perturbation parameter? This question can be answered by considering a 
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numerical example based on the solution presented in equations (76). Consider an alternate 
choice for <5C in which the number of repetitions of a 1 has been reduced by one. A solution for 
the uncertainty dynamics is 


<5C = diag a 1 , a' , t/, c' 



10 10 
1 s 0 1 


(77a) 


(s) 


1 


■s' 2 + ds + e 


0 0 0 

1 1 1 


L e< , (s) — [0] 4x 4 


(77b) 


with an appropriately reduced sector radius matrix R. 

One would expect the uncertainty dynamics in equations (76) to result in a more conservative 
uncertainty bound than that from equations (77) since, in equations (76), a' can effectively 
take on three different values simultaneously rather than only two, as in equations (77). The 
conservatism associated with the added degree of freedom is demonstrated in figure 9. These 
bounds are generated with equation (58) 3 , with the values for the system parameters given by 


a = 5.4 b = 10.0 c = 25.25 d = 10.0 e = 100.0 
r a = 4.5 r b = 2.0 r c = 5.0 

Even though the bound obtained from equations (77) is less conservative than that from 
equations (76), it still exhibits conservatism. Some of the conservatism is because a 1 can still 
simultaneously take on two different values inside its sector, when in reality it can only have one 
value. Additional conservatism is introduced when a' is actually real valued, since the sector 
bound assumes complex perturbations. This conservatism is evident when the sector bounds 
are compared with an exact bound for the assumption of only real-valued perturbations, as in 
figure 9. 


Figure 9. Conservatism of uncertainty bound due to perturbation representation. 


3 A listing of the sector bound algoritlm is presentedin appendk B. 


19 



The exact bound is obtained by using a constrained optimization procedure (Jacob 1972). 
The procedure uses a nongradient search method to maximize the singular value of the additive 
uncertainty matrix, A in figure 5, associated with the uncertain parameters. In this case A 
corresponds to the transfer function matrix between d(s) and u(s) (eq. (63)). The uncertainties 
are represented as real perturbations about the nominal values of the uncertain parameters. The 
maximum variations are symmetric about the nominal values and equal to the sector radii. 

As shown above, conservatism is introduced into the sector bound when a perturbation 
parameter is repeated more times than required. Therefore, it would be desirable to be able 
to identify when a parameter is excessively repeated so that the conservatism can be reduced. 
Alternatively, it would be useful to have a method by which the correct number of parameter 
repetitions could be determined before the uncertainty dynamics are obtained. In that way 
conservatism due to excessive repetitions could be eliminated. 

In the above example, the solutions of equations (71a) clearly show that the third appearance 
of a' does not affect the system response. Recall, however, that the solutions are not unique. 
Not all solutions so readily show that a' is excessively repeated. For example, if the “inspired 
guesses” in equations (74) are chosen so that 


L^.;, (s) 


1 

s' 2 + ds + e 


1 

0 


0 10 0 

10 11 


then one solution for L e . u (s) is, 



k i 0 1 - k i 1 0 

k -2 s 1 — k -2 0 1 


(78) 


(79) 


where k\ and &2 are arbitrary polynomials in s. 

If k\ — k 2 — 1, the third appearance of the perturbation a 1 is not excited by any system 
input and, therefore, has no effect on the system response. Similarly, if k\ — k 2 — 0, the first 
appearance of a' has no effect on the system response. For other values of k\ and k 2 , the fact 
that a 1 only needs to appear twice is not obvious since each appearance of a' is excited by the 
system input u(s) and generates a perturbation to the system output d(s). Therefore, it is not 
appropriate to use zero rows or columns alone to identify excessively repeated perturbations. 

At present there is no reliable test using properties of the L^s) and L eUi (s) matrices to 
identify excessively repeated perturbation parameters. The recommended approach is to start 
with a single independent parameter and sequentially add repeated parameters until the system 
of equations (e.g., eqs. (71)) can be solved, subject to the restrictions that the solutions be 
minimal polynomials. 


5.3 Denominator Uncertainty for a Multivariable System 

The approach to solving the denominator uncertainty problem is identical to the approach 
used for numerator uncertainties. The differences occur in the specific equations that need to 
be solved. The next example addresses generic denominator uncertainties to highlight these 
differences and to address the properties of the resulting system representation. 

The system in equation (60) is considered once again, but in this case each of the denominator 
coefficients are uncertain and known to lie in sectors. That is, 

d^d + d' e^e + e' (80) 
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where d and e are the nominal values of d and e and d' and e' are the perturbations of d and 
e. The sector bounds on the coefficients are such that the perturbation is centered at the origin 
with radius r (e.g., d' is in a sector centered at zero with radius ) . 

In this case the system frequency-domain input-output equation can be written as 


(s 2 + ds + e)y (s) 


s + a a 

u (s) - 

d's + e' 

0 1 

s + b s 2 + as + c 

0 

d's + e' 


The perturbation response is, therefore. 


d(s) 


A 


£ (£) 

*0) 


-1 

s 2 + ds + e 


d's + e' 

0 


0 

d's + e' 


y 0) 


(81) 


(82) 


where <5(s) is the denominator polynomial and r(s) is the matrix of numerator polynomials. 
From the block diagram in figure 4, equation (23a), and equation (82), one finds that two 
equations must be satisfied by the choice of uncertainty dynamics matrices. These equations are 


where 


T (s) N y U (s) _ 

S (s) d yu (s) 

Nyv ( s ) fig Neu (s) 
dyv d eu ( 5 ) 

(83a) 

T (s) N y V (s) _ 
6 (s) d yv (s) 

N yv ( 5 ) ^ N ev (5) 
dyv (^) d^v (s) 

(83b) 

T / w,'. _ ( S ) 

"" ) dy U (s) 

( ] dy V ( S ) 

(84a) 

N ett (s) 
'“*■ d eu ( s ) 

N et ,(s) 
" ( )_ d„.(s) 

(84b) 


and where the N(s) terms represent matrices of numerator polynomials and the d(s) terms 
represent denominator polynomials. 


Choose d yv (s) — S(s) — d yu (s), consistent with the previous example, and choose d ev (s) — 
d eu (s) — S(s). In this case equations (83) simplify to 


-T(s) ! > I - S ,/,■ ! > I SC N eu (s) (85a) 

-r (s) N yv (s) - Ny V (s) 6C N ev (s) (85b) 

As in the previous example, the solutions of equations (85) are not unique and, in order 
to make the solutions more tractable, the same restrictions are placed on the elements of the 
matrices N 2/ . u (s), ’N yv (s) , N ett (s), and N et ,(s). Thus, the elements of N(s) are required to be 
minimal polynomials (i.e., they only have terms involving either nonnegative integer powers of 
s or zero and are of the lowest possible order and the least possible number of terms). 

The first step in solving equations (85) is to choose <5C and its associated sector radius 
matrix R. In this example, each uncertain parameter in <5C must be repeated once, as if it is 
associated with two independent sectors. The repeated perturbation parameters d' and e' are 
required to characterize the overall effects of the uncertainty. Therefore, for this example, 

SC - diag[d', d', e', e'] R = diag [r d , r d , r e , r e ] (86) 


21 



One should expect the resulting sector bound to have some conservatism associated with the 
repeated perturbation parameters. 

Next consider the uncertainty dynamics numerator matrices. For dimensional consistency 
N y V (s) must be a 2 X 4 matrix, N e . u (s) must be a 4 X 2 matrix, and N et ,(s) must be a 4 X 4 
matrix. Represent the numerator matrices as 


Nyu (s) 


s + a a 


k n 

l 

CM 

s + b s 2 + as + c 


>’21 

&22 J 


(87a) 


r (s) 

... r n r 12 
_r 2 i r 22 . 

2 

CS 

11 

k 1 

^11 

^21 

^12 

^22 

^13 

^23 

— 1 

to I— 1 

1 1 






" **11 

**12 

**13 

**14 

m n 

m 2 i m .31 ***41 

Nf ii (s) 

— 

**■21 

**22 

**23 

**24 

***12 

***22 in-32 ***42 . 


**31 

**32 

**33 

**34 





- **41 

**42 

**43 

**44 


(87b) 


(87c) 


Substituting equations (86) and (87) into equations (85), expanding the elements, and 
grouping like terms results in a system of equations for each uncertain parameter. For this 
example, the two systems that result can be written in an indicia! form as follows. 


System for d: 


System for e: 


skij — hl m lj T h'2 m 2j 
slij — ^'l**l 2 T ^2**2j 


{i,j= 1,2)1 
(*' = 1,2; j = 1, 2,3,4) j 


( 88 ) 


hi] — h3 m 3j + ^'4***4j 

hj — ^’3**3 j T ^'4**4j 


{i,j= 1,2)1 
(* = 1,2; 1.2.5. 1) J 


(89) 


Again, inspired guesses are used to solve these two systems of nonlinear equations. Let 


^11 — ^22 — 1 h2 — /21 — 0 ^13 — ^24 — 1 ^14 — ^23 — 0 (90) 


The remaining numerator elements are determined next, subject to the imposed limitations 
discussed above. As a result the: complete uncertainty dynamics representation becomes 


<5C = diag^, d' ,e' , e ; ] 


(91a) 


L^O) 


s“ + ds + e 


s(s + a) s(s + b) 


s + a 


sa 


s(s“ + as + c) 


s + b 

s 2 + as + c 


L ^’ (S) - s 2 + ds+e 


Lec(s) 


1 


s 2 + ds + e 


10 10 ' 

0 1 0 1 

s 0 s 0 ' 

0 s 0 s 

10 10 
0 10 1. 


with the associated sector radius matrix R from equations (86). 


(91b) 

(91c) 


(91d) 
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An example of the sector bound for the denominator uncertainty case (eq. (58)) is shown 
in figure 10. The same set of nominal values for the system parameters used in the previous 
example are used here as well; however, in this case the uncertain parameters are d and e. The 
nominal values and the radii for the uncertain parameter sectors are given by 

a = 5.4 b = 10.0 c = 25.25 d=10.0 

e= 100.0 r d = 5.0 r e = 40.0 

Notice the conservatism associated with the need to treat the two uncertain parameters 
as four independent sectors when compared with the “exact” bound (based on real-valued 
perturbations). The exact bound is computed with constrained optimization, as in the previous 
example. 

It is also interesting to note that one of the sufficiency conditions for existence of the sector 
bound (namely, <j[RL et ,(s)] < 1) is violated for values of r e greater than 50. In order to compute 
a sector bound for r e > 50, one could scale the inputs and outputs of the uncertainty dynamics 
6L 2/ . u (s) to reduce the norm of R. 


Figure 10. Sector-based uncertainty bound for denominator uncertainty. 

5.4 Multiple Parameter Uncertainties in a Second-Order System 

The solution of the sector bound for a general system is accomplished by combining the 
solution methods for numerator and denominator uncertainties. A general representation of the 
uncertainty disturbance d(s) can be represented as a function of the inputs and outputs of thq 
true system: 

d (s) = Ti (s ) u (s) - T 2 (s ) y (s) (92 ) 

The transfer matrices relate the inputs and outputs of the true system to the uncertainty 
disturbance. The uncertainty disturbance expression in equation (92) can be combined with 
equation (23a) to obtain a relation in terms of u(s) and v(s): 

d (s) = [r 1 (s) - r 2 (s) L yu (s)] U (s) - r 2 (s) L^, (s) v (s) (93 ) 
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The uncertainty disturbance can be written in an alternate form by inspection of the block 
diagram depicted in figure 4: 

d (s) - L yv (s) 6C [L e „ (s) u (s) + L e; , (s) v (s)] (94) 

Equating the two expressions for the uncertainty disturbance (eqs. (93) and (94)) and grouping 
the u-terms and v-terms results in two expressions that the uncertainty dynamics must satisfy: 


r i (s) r , (s) L yu (s) — Ly V (s) 6C L e „ (s) 


(95a) 


-r 2 (s) Lyp (s) = L yp (s) 6C L e< , (s) (95b ) 

Note that the two expressions in equations (95) are algebraic in nature. They correspond 
to a set of scalar equations that are nonlinear in the unknown transfer function elements. The 
solutions of these equations are similar to those of the special cases discussed previously. The 
solutions are not unique because the equations underspecify the solution. Analogous to the 
previous examples, however, a solution is readily obtained by imposing the same practical 
restrictions on allowable forms of the solution. 

In order to see how the generalized solution approach can be applied to a problem with 
physical significance, consider the transfer function in equation (96). The nominal system is 
unity but has unmodeled dynamics that are known to be in the form of a second-order damped 
oscillator. 4 This system has a true transfer function of the form 


y 0) 

u (s) 


G(s) 


IcW 
I O/V, < 


— Gq (s) + A (s) — 1 — 


' + 2 


(96) 


where the frequency ui and damping ratio £ are uncertain but are bounded as follows: 


0.4 < ( < 0.8 


22 < u < 27 rad/sec 

A sector bound for this example is computed two ways. The first way might be termed the 
naive approach, in which each appearance of an uncertain parameter is treated as an independent 
uncertainty. The second solution approach makes use of the physical structure of the system 
to recognize that the same parameter appears in both the numerator and denominator of the 
system transfer function. One would expect the second approach to give a less conservative 
representation of the uncertainty, and this will be shown to be the case. 

For the first case the system is modeled as shown in equation (97): 


Gi(s) 


a 

s 2 + bs + c 


(97) 


The parameters a, b, and c are independent and related to u> and £ in such a way that the 
nominal parameter values and sector radii are 


a = 606.5 b = 30.4 c= 606.5 

r a = 122.5 r b = 12.8 r c = 122.5 


That is, a and c lie in sectors centered at 606.5 with radii 122.5, and b lies in a sector centered 
at 30.4 with radius 12.8. 

4 This example is adapted from the examplein Stein and Doyle (1978). 
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The perturbation matrix <5C and the T matrices are given by 


<5C = 


0 

0 


0 0 
b' 0 
0 d 


ri(s) 


■s' 2 + bs + c 


r 2 («) 


b's + c' 
s 2 + bs + c 


(98) 


where the diagonal structure of the <5C matrix is dictated by sector theory. One solution 
of equations (95) for these values of <5C, T^, and T 2 , with Gi(s) substituted for L 2/ . u (s), is 
summarized in equations (99): 


L</(.(s) 


1 


■s' 2 + bs + c 


[ 1 , 1 , 1 ] 


(99a) 



1 


s 2 + bs + c L 


T bs T c, — 3s, — ci] 


(99b) 


L e ;i(s) 


1 

s 2 + bs + c 



0 0 

— s — s 

-1 -1 


The sector radius matrix for this case is 


R = 


r a 

0 

0 


0 

rb 

0 


0 

0 

r c 


(99c) 


( 100 ) 


Notice that L^s), L et ,(s), and L eUi (s) have the same denominator corresponding to the 
characteristic polynomial of the nominal plant. This is not required but significantly simplifies 
the solution and the structure of the uncertainty dynamics. This particular solution, once the 
denominators are chosen, is determined by the choice of the elements of the numerator of the 
Ly(,(s) matrix. The remaining numerator terms are then specified by the governing relations, 
equations (95). 

Also notice that the upper row of the L e; ,(s) matrix, the one associated with the uncertain 
numerator parameter a, is zero. This is because the L et ,(s) matrix is the feedback block inside 
the uncertainty dynamics block diagram (fig. 4) and is required only when there is uncertainty 
in the denominator parameters. Since a only appears in the numerator of Gi(s) no feedback 
around the <5C matrix is needed. 

The sector-based uncertainty bound corresponding to equation (54) for the above solution 
(three independent parameters) is presented in figure 11 . The “absolute” uncertainty bound 
determined by constrained optimization methods is also presented for comparison. Notice the 
significant conservatism associated with the sector-based bound, especially at low frequencies. 
This conservatism can be attributed to several factors, as discussed previously. The source of 
most of the conservatism, however, is the treatment of the parameters as three independent 
uncertainties. This is demonstrated in the next example. 


25 



Figure 11. Comparison of uncertainty bounds. 


Consider the system modeled as shown in equation (101), in which the parameter a appears 
in both the numerator and denominator. This accounts for the fact that in the actual system,, 
equation (96), this is indeed the case. The nominal parameter values and sector radii for a and 
b are the same as for the previous case. 


G 2 («) = 


s 2 + bs + a 


( 101 ) 


The perturbation matrix <5C and the T matrices for this case are given by 

b's + a' 


6C = 


a' 0 
0 b' 


Ti(s) 


s 2 + bs + a 


r 2 («) 


■s' 2 + bs + a 


( 102 ) 


A solution for the uncertainty dynamics for the two-parameter case, shown in equations (103), 
is very similar to the three-parameter solution, with the exception that equations (103) are of 
lower order: 


L^„ (s) 


1 


■s' 2 + bs + a 


[ 1 , 1 ] 


(103a) 



1 

s 2 + bs + a 


s - + bs, — as 


(103b) 


L e ,, (s) 


1 

s 2 + bs + a 




The sector radius matrix for this case is also of lower order: 


(103c) 


R = 


r a 

0 


0 

rb 


(104) 
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The sector-based uncertainty bound corresponding to equation (58) for the case of two 
independent uncertain parameters is presented in figure 11 along with the “absolute” uncertainty 
bound and the three-parameter bound. Notice that this bound is considerably less conservative 
than that of the previous case. This is attributable to the fact that the physical structure of the 
true system is more accurately represented. Some conservatism remains, however, because of 
the reasons cited previously. Conservatism is also introduced because the parameters a and b, 
while treated as such, are not truly independent of one another, since both depend on the value 
of ui. 

The examples presented above were used to demonstrate the manner in which conservatism 
enters the sector-based uncertainty bound. In addition, the examples led the development of 
a general form for the equations that govern the uncertainty dynamics 6 L yu (eqs. (95)). The 
nonuniqueness of the solution to the governing equations is one source of conservatism that was 
demonstrated in the first example. The need to repeat single uncertain parameters multiple 
times is another source of conservatism that was demonstrated in the second example. The need 
to treat multiple dependent parameters as independent is yet another source of conservatism 
that was demonstrated in the third example. The need to repeat some parameters and treat 
dependent parameters as independent can in most instances be traced to the requirement that 
the perturbation matrix <5C be diagonal. 

The examples have also shown that the sector-based uncertainty bound provides a means 
to readily compute^jin closed form, a frequency-dependent characterization of the influence of 
individual uncertain parameters on the overall frequency response of an LTI-MIMO system. 
This characterization of the system uncertainty can be directly applied to robustness tests and 
robust control design tools. 

6 Concluding Remarks 

This paper has addressed several aspects of applying sector stability theory to the problem 
of characterizing plant model uncertainty. Sector properties were applied to linear multivariable 
systems with Safonov’s approach. A general matrix representation for linear, time-invariant, 
multivariable systems was presented to facilitate the application of the sector properties. The 
resulting sector-based uncertainty bound was combined with a representative singular-value 
robustness condition to obtain a useful robustness measure that accounts for known uncertainty 
in the plant model. Finally, the sector-based approach was applied to several forms of parameter 
uncertainty to highlight the specific issues involved in applying the sector-based approach to 
linear, time-invariant, multivariable systems and to identify properties of the sector bound 
associated with fundamental forms of plant model uncertainty. 

Based on the results obtained, the characteristic advantages and limitations of the sector- 
based approach are clarified. Three primary advantages of the sector-based approach become 
evident: the approach provides a way to effectively describe overall plant model uncertainty when 
bounds on individual model parameters are known, the approach is multivariable by nature and 
so directly applies to linear multivariable dynamical systems, and the approach is compatible 
with established singular-value-based robustness measures. 

The major limitations associated with this approach are the conservatism of the resulting un- 
certainty bound and the difficulty in solving the systems of nonlinear equations that determine 
the uncertainty dynamics. Conservatism may be introduced by the diagonal matrix representa- 
tion of the uncertain parameters required to apply sector concepts to the uncertainty modeling 
problem. Each uncertain parameter is assumed to be independent of all other uncertain param- 
eters. Conservatism is introduced whenever repeated parameters are required. Mathematical 
conservatism is introduced by the process of merging the sector bound with singular-value con- 
cepts. This conservatism can be attributed to two main factors: the merging process requires 
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the use of inequalities, and the uncertainties are assumed to be complex valued with an arbitrary 
phase angle (i.e. , bounded only in magnitude). Another factor that may have some effect on 
the degree to which the sector bound is conservative is associated with the nonuniqueness of the 
uncertainty dynamics representation. The systems of nonlinear equations associated with the 
uncertain parameters can be readily solved, as shown in the examples, but they are underspeci- 
fied. The level of conservatism is increased if the solution of the nonlinear equations results in 
excessively repeated uncertain parameters. 

A feature of the methodology that has not yet been exploited effectively is the nonuniqueness 
of the solution to the uncertainty dynamics equations. There may be a “best” solution that 
minimizes the overall effect of the various forms of conservatism on the accuracy of the sector- 
based uncertainty bound. The additional degrees of freedom might be able to be used to generate 
such a solution, which results in the least conservative sector bound. Solution methods that 
exploit the sector approach and result in the least conservative uncertainty bound are subjects 
for future research. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
January 21, 1992 
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Appendix A 
Proofs 

Proof of Lemma A1 

Lemma Al. Given two complex matrices A and B such that A is n X n, A* = A, and 
A > 0 (i.e., positive definite), B is r X n, rank (B) = r, and r < n, 


then, 



BA _1 B* 


A > B*AB 


(i.e., x*Ax > x*B*ABx). 


(Al) 

(A2) 


Proof. Define 


Bt ^ A _1 B* 


BA _1 B* 


-l 


which is the right pseudoinverse of B with respect to A. Then, 


(A3) 


A = 


BA _1 B* 


i-l 


BA _1 B* 


BA _1 B* 


i -l 


BA _1 B* 


BA _1 B* 


i -l 


-l 


BA _1 AA _1 B* 


BA _1 B* 


i -l 


B 1 ] AB f 


(A4) 


Now consider 


A - B*AB 


x = x*Ax - x*B*ABx 

* A - 


x*B* 

B f 

*AB f Bx 

X* 

B f B 

* 

A 

B f B 


(A5) 


where x is an appropriately dimensioned complex vector. In addition, since BB^ = I (by 
definition) and since 


B f B 

= A -1 B* 

BA _1 B* 

1 BA _1 B* 

BA _1 B* 

-l 

B 


= A _1 B* 

BA _1 B*1 

-1 

B 




= B f B 


then B^B is a projection such that 


B f B 


X = X 


J ? _L (B) 


(A6) 

(A7) 


where is a vector that lies in the orthogonal complement of the null space of B. Both the 

null space of B and its orthogonal complement are spaces of lower dimension than those of A. 
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Any vector x can be written as the sum of two components, one that lies in the null space of 
B, 17(B), and one that lies in its orthogonal complement, so that 

x = x„ (B ) + x,± (B) (A 8 ) 

As a result, 


x* (A - B*AB J x = x*Ax - Xj?±(B) A Xj?±(b) 

= X ^(B) Ax »?(B) (A9) 

Since A is given to be positive definite and since x^ B ^ lies in a space of lower dimension than 
A, then 

A — B*AB > 0 (A10) 


(positive semidehnite) and so 

A > B*AB 


(AH) 


Sector Bound Existence Tests 

In this section simple tests to check for the existence of a system sector bound are developed 
and proven. In order to calculate a system sector bound several matrices need to be formed. 
Some of these matrices require the inverses of other matrices to exist. It is advisable to check 
for the existence of these matrix inverses before attempting to determine the sector bound itself. 
Two tests are developed to determine the existence of these matrix inverses. 

Existence test for Q” 1 . The inverse of Q„ must exist in order to calculate the matrices in 
equations (40). For a matrix to have an inverse it must be definite. Since Q„ is Hermitian it 
must be positive definite. This fact is used to develop the existence test. 

Since Q„ must be positive definite, 


x*Q„x > 0 (A12) 

O x* [I- L*„R*RL e „]x > 0 (A13) 

O x*Ix - x*L* 0 R*RL eo x > 0 (A14) 

Now if x*x 7 ! 0, then 

x*LT,R*RL e „x 

O ^ — < 1 A15 

x*x 

Using the Euclidean norm notation gives 


O 


|RL P 


|2 
I E 


|2 
I E 


< 1 


(A16) 


which is guaranteed to be true if the following singular-value inequality is satisfied: 


u (RL e? j) < 1 

Therefore, equation (A17) is a simple test for the existence of Qif 1 . 


(A17) 
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Existence test for Q ra . Since Q ra is defined to be the inverse of a matrix product, that 
inverse must exist. Recall that 

Qra — (A18) 

For Q ra to exist, Q” 1 must exist and L^Q^L*,, must be positive dehnite. Assume that the 
existence of Q” 1 has already been established. 

Since Q” 1 exists, it must be full rank and it can be factored (e.g., by Cholesky decomposition). 
Let the factors be defined by 

Q- 1 = FF* (A19) 

Therefore, it must be shown that 

\j yv FF*\jy V > 0 (A20) 

Note that Q„ is square and let its dimension (and rank) be n. Then the dimension (and rank) 
of F is also n. Note also that L yv is not necessarily square; let its dimensions be to X n. Clearly 
then, L^FF*!^ is an to X to matrix. Such a matrix is only positive dehnite if (1) m < n and 
(2) rank(L 2/ „) = to. 

If to > n, then L^FF*!^ is deficient because the rank of L yv can be at most n and so the 
rank of L^FF*!^ can be at most n as well. Also, since the matrix product L^FF*!^ is 
Hermitian it will be positive dehnite if it is full rank. Therefore, if condition (1) is satished and 
the rank of L yv is to, then L^FF^*,, is of full rank. 

Both of the tests presented above can be readily performed numerically. They should be 
performed as part of the algorithm used to determine system sector bounds to ensure that the 
problem is properly formulated and that the bound resulting from such an algorithm is valid. 
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Appendix B 

Sector Bound Algorithm 

The sector algorithm is coded with command hies for MATRIX^, 5 a commercially available 
linear systems analysis and design package. Three separate command hies are used: safsys.dat 
is used to set up the uncertainty dynamics matrices, safbnd.dat is used to set up the sector 
bound matrix, and calcs.dat is used to actually compute the sector-based uncertainty bound. 
An additional command hie, logdist.dat, is used to create a frequency vector with points evenly 
distributed on a logarithmic scale. 

safsys.dat 

The command hie safsys.dat is problem specific and needs to be customized for each solution 
of the uncertainty dynamics matrices (i.e., L^, L^,, L e(I , and L eu ). A listing of safsys.dat for 
the denominator uncertainty example from section 5.3 is presented below. 


It 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

a 


Command File - safsys.dat 


Author 

Martin R. Waszak 

Date 

6/29/90 

Description : 

Formulate the representation of the system for use 
in calculating Safonov’s sector uncertainty bound. 


The formulation is problem dependent and must be reprogrammed 
for every new problem to be addressed. 


This formulation only applies for the denominator uncertainty 
example case which appears in section 5.3. 

5.4; b = 10.0; c 

= 25.25; d = 10.0; e = 100.0; 


svar = jay*omega(i); 

den = svar**2 + d*svar + e; 


ns = 2;nv = 4; 

gO = 1/den*[svar+a, a; svar+b, svar**2+a*svar+c]; 

gyu = 1/den*[svar+a, a; svar+b, svar**2+a*svar+c]; 

geu = 1/den*[ svar*(svar+a), svar* (svar+b), svar+a, svar+b;... 

svar*a, svar*(svar**2+a*svar+c), a, svar**2+a*svar+c]’; 
gyv = 1/den*[ 1 0 1 0; 0 1 0 1 ]; 

gev = 1/den*[ svar, 0, svar, 0 ; 0, svar, 0, svar;1, 0, 1, 0; 0, 1, 0, 1]; 

clear svar den a b d c e 
return 


® MATRIX ^ is a registered tradoiame of IntegratedSystems Inc. 
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safbnd.dat 


The command file safbnd.dat is problem specific and needs to be customized for each choice 
of the perturbation matrix <5C. The sector radius matrix is determined from the individual 
parameter sectors associated with the uncertain parameters represented in <5C. A listing of 
safbnd.dat for the denominator uncertainty example from section 5.3 is presented below. 


// 

// 

// 

It 

II Author 
// Date 

// 

// Description 

// 

// 

// 

// 

// 


Command File - safbnd.dat 


Martin R. Waszak 
6/29/90 

Formulate the sector radius matrix associated 

with the sector bounds on each of the uncertain elements. 

This formulation only applies for the denominator uncertainty 
example case in section 5.3. 


r = diag([5. 0,5.0, 40. 0,40.0]); 
return 


calcs.dat 

The command hie calcs.dat computes the sector-based uncertainty bound described in 
equation (58). It applies for all forms of uncertainty dynamics (e.g., additive, multiplicative 
reflected at the plant input, and multiplicative reflected at the plant output). A listing of 
calcs.dat is presented below. 


// 

// 

// 

// 

// Author 
// Date 

// 

// Description 

// 

// 

// Remarks 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 


Command File - calcs.dat 


Martin R. Waszak 
6/29/90 

Calculate a singular value bound on model 
uncertainty based on sector stability theory. 

This algorithm can determine sector bounds, for both 
modeled parameter uncertainties and for unmodeled dynamics 
with known structure but with uncertain parameters, where 
the uncertainties can be represented by sectors. 

The variable names do not correspond to the matrices 
defined in the theoretical sector formulations. 

They do have some relationship with the ‘standard’ 
notation, however. They have been altered 
to improve the computational aspects of the algorithm. 
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// 

Inputs 

logdist.cmd 

- command file to formulate frequency vector 

// 



evenly spaced on a log scale 

// 


omega 

- vector of frequencies (from logdist.cmd) 

// 

// 

// 


n 

- number of frequencies 


safsys.dat 

- command file to formulate 

// 



system representation 

// 


gO, gyu, gyv, 


// 


geu, gev 

- system matrices (from safsys.dat) 

// 


ns 

- size of the nominal system (from safsys.dat) 

// 


nv 

- number of component uncertainties 

// 

// 

// 



(from safsys.dat) 


safbnd.dat 

- command file to formulate sector 

// 



component uncertainty bound matrix 

// 


r 

- matrix of sector bounds 

// 

// 

// 

// 

// 



(from safbnd.dat) 

Outputs 

bs 

- system sector bound 

Modified 

7/7/88 

- to correct the expressions for ‘tl’ and 12’ 

// 


4/1 8/89 

- to modify the sector bound equation; bs(i) 

// 


7/1 2/89 

- to calculate svd(r*gev) for existence test 

// 


6/29/90 

- to apply for additive uncertainties 


//======================================================== 

// 

// Obtain frequency vector evenly spaced on a log scale 

// 

wmin=0.1 ;wmax=1 00;npts=300;exec(‘logdist.cmd’) 
n=300; 

// 

// Perform the bound calculation for each frequency (FOR LOOP) 

// 

display(‘Forming system matrices and bounds - (this will take a while)’),... 
for i=1 :n;... 

...II 

...II Obtain the special sector system representation 
...II and matrix of element sector bounds 
exec(’safbnd’); exec(’safsys’);... 

...II 

...II Calculate the necessary matrices 
tl = r*gev; t2=r*geu;... 

...II Calculate svd(r*gev) which should be < 1 for existence 
sv1=svd(t1);... 

if svl (1 ) > = 1.0 then display(‘existence test violated’), return, end,... 
ml = inv( eye(nv) - tl'*t1 ); m2 = inv( eye(nv) - t1*t1 ');... 
q =real( gyv*m1*gyv'); p = real( t2'*m2*t2 );... 
gc = gyu + gyv*m1*t1’*t2;... 

...II 

...II Calculate Cholesky factors 
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cfq = chol(q); cfp = chol(p);... 

...II 

...II Calculate singular values of appropriate matrices 

svl = svd(gc - gO);... 

sv2 = svd(cfq); sv3 = svd(cfp);... 

...II 

...II Calculate sector bound and then loop back for another frequency 
bs(i) = ( svl(l) + sv2(1)*sv3(1) ) ; end 

// 

display('AII done! - sector bound stored in ”bs” ’) 

// 

clear gO gc gyu gyv geu gev ns nv ml m2 tl t2 n p q r 

clear cfq cfp svl sv2 sv3 sv4 

return 


logdist.dat 

The command file logdist.dat is used to generate a logarithmically distributed frequency vector 
in a range between specified minimum and maximum values. 


// 

// 

// 

// 

// Author 
// Date 

// 

// Description 

// 

// 

// 

// 

// Inputs 

// 

// 

// 

// 

// 

// 

// 

// Outputs 

// 

// 

// 

// 

// Remarks 

// 

// 

// 

// 


Command File - logdist.cmd 


Martin R. Waszak 
4/12/88 

This function calculates a range of points that 

are uniformly distributed on a loglO scale. It is ideally 

suited for use in plotting things on a log scale 

(for example - bode plots, singular value plots, etc.). 

The following temporary variables must be defined 
in your matrix-x stack. They will be cleared from the stack 
after the vector of points is formed. 

npts = # of data points 

wmin = minimum value of the range of values 

wmax = maximum value of the range of values 

The range of logarithmically distributed 
values will be output in the following vector - 

omega = vector of range of values 

Since this is a command file, the variables used 
here are global (not local) and care must be taken that 
these variable names are unassigned in your calling 
program. 
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6/29/88 


// Modified 

// 

// 

// 

// 

//========= 

// 


- to display notification that calculation has 
begun 

6/30/88 - to improve documentation and improve 
generality 


display(‘Calculating a vector of points evenly spaced on a log scale.’) 

// 

for i=1:npts; omega(i) = wmax**(i/npts) * wmin**((npts-i)/npts); 

// 


display(‘The logarithmically distributed range of -’) 
npts 

display(‘values with minimum value -’) 
wmin 

display(‘and maximum value -’) 
wmax 

display(‘is stored in the column vector “omega.” ’) 

display('Note : wmin, wmax and npts have been cleared from the stack.’) 

clear wmin wmax npts 

return 
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Nomenclature 

A arbitrary matrix (complex or real valued depending on application) 

B arbitrary matrix (complex or real valued depending on application) 

C(s) sector center matrix 

c sector center; uncertain parameter 

d uncertain parameter 

d(s) perturbation response 

d eu (s) denominator polynomial of L eUi (s) 

d ev (s) denominator polynomial of L et ,(s) 

dy V (s) denominator polynomial of (s) 

e uncertain parameter 

e(s) input to perturbation matrix <5C 

F(s) matrix of complex- valued functions 

/(•) function of (•) 

G(s) nominal plant transfer matrix 

K(s) feedback compensation 

fcj arbitrary polynomial 

kij i,j element of N^ ((s) 

Leu(s) transfer matrix that transforms the vector of system inputs 

L e .;,(s) transfer matrix that allows parameters in <5C to appear in denominator of the 

uncertainty dynamics 

L 2 / . u (s) transfer matrix for the nominal plant model 

L^(,(s) transfer matrix that transforms the output of the perturbation block into the 

perturbation response 

lij i, j element of N yv {s) 

m-ij i, j element of N eUi (s) 

N e „(s) numerator of L eUi (s) 

N y V (s) numerator of L y ( s ) 

n-ij i, j element of N et ,(s) 

P„ defined in equation (44b) 

P ( , defined in equation (34) 

Q„, defined in equation (44c) 

Q ( , defined in equation (28a) 

R(s) sector radius matrix 

R ( , defined in equation (28b) 

r sector radius 
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r (») 

s v 

S 

T(s) 


t 

u («) 

u(t) 

v(s) 

X 

y(«) 

A 

A 

(5(s) 

<5C 

6c 

6hy U 

r H (s) 

a(«) 

A(«) 

U) 

c 

(•) 

(•) 

(•)' 

(•)* 

II • lb 
II • We 
<(•).(•)) 


sector radius for parameter (•) 
defined in equation (28c) 

Laplace variable 

linear, time-invariant, multivariable operator 

defined in equation (44a) 

defined in equation (30) 

independent variable 

plant input vector 

arbitrary real function of t 

output from perturbation matrix <5C 

arbitrary vector compatible with A 

plant output vector 

uncertainty representation 

second term on right side of equation (44a) 

uncertainty dynamics denominator polynomial 

perturbation matrix 

perturbation from sector center 

uncertainty dynamics 

matrices of uncertainty dynamics numerator polynomial 
maximum singular value of (•) 
minimum singular value of (•) 
frequency 

uncertain parameter; damping ratio 

Fourier transform of (•) 
nominal value of (•) 
perturbation of (•) 
complex conjugate transpose of (•) 
absolute value of (•) 

2-norm 

Euclidean norm of (•) 
inner product 
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